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Introduction 

At sufficiently high density, the transition to deconfined strange quark matter is widely 
expected. Chiral symmetry is spontaneously broken in the QCD vacuum. Lattice QCD 
simulations at nonzero temperature T and zero baryon-chemical potential jiB indicate that 
chiral symmetry is restored above a temperature T ~ 150 MeV [1]. The NJL model is 
an effective theory which is believed to be related to QCD at low energies, when one has 
integrated out the gluon fields. The NJL model might yield reasonable results at the density 
range where confinement is no longer crucial but chiral symmetry as a symmetry of full QCD 
remains to be important. The NJL model has proved to be very successful in the description 
of the spontaneous breakdown of chiral symmetry exhibited by the true (nonperturbative) 
QCD vacuum. This model has been extensively used over the past years not only to describe 
hadron properties [ ] (see for reviews [4, 5]) and phase transitions in dense matter [6, 7, 8] 
but to describe the quark strange stars as well [9, 10, 11, 12, 13]. The detailed properties of 
the quark phase in compact stars has been a topic of recent interest [14] (for a review see 
[15, 16]). 

Quark strange stars are astrophysical compact objects which are entirely made of 
deconfined u,d,s quark matter {strange matter) staying in /3-equilibrium. The possible 
existence of strange stars is a direct consequence of the conjecture [17] that strange matter 
may be the absolute ground state of strongly interacting matter. 

The three-flavor NJL model has been discussed by many authors, e.g. [ ! 8] . For the quark 
phase we follow Buballa and Oertel [ ] in using the three-flavor version of the NJL model. 

The aim of this paper is to investigate a strange quark star within the enlarged SU(3) 
Nambu-Jona-Lasinio (NJL) model. The comparison of a quark star properties obtained in 
the Quark Mean Field (QMF) approach to a neutron star model constructed in framework the 
Relativistic Mean Field theory (RMF) will be made. This paper is organized as follows. 
In Section 1 there are presented general properties of the NJL model in the mean field 
approach based on the Feynman - Bogolubov inequality for free energy of the system. In 
this section the employed equation of state (EoS) is calculated for NJL model. The EoS is 
used then to determine the equilibrium configurations of the quark star in Section 2. Finally, 
in Section 3 the main implications of the results are summarized. 

Namby-Jona-Lasinio model 

The NJL[2] model was widely used for describing hadron properties ["] and the chiral 
phase transition [ i ]. The enlarged (ENJL)[ ] simplest version of the model is given by 
the Lagrangian: 

L = q{zYd, - mo)q + ^G. j^HqX^'q)' + (gA"^75g)'] - [] iWf) 

a=0 f={u,d,s} 

-Ig.T. [{qY^'^q) {qi^X^^q) + (g7^5A"g) (g7M75 A"g)] (1) 

^ a=0 
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+ ^ E 'L'n'd^Lf - mfLfLf + Bo 
/=i /=i 

The first term contains tlie free kinetic part, including the current quark qf = {qu,Qd,Qs} 
masses niQ which break explicitly the chiral symmetry of the Lagrangian and the term 
representing the free relativistic leptons Lf = {e~, The fermion fields are composed of 
quarks and leptons (electrons, muons). Here q denotes a quark field with three flavors, u, d 
and s, and three colors.We restrict ourselves to the isospin SU(2) unbroken symmetric case. 



rrin 



mi, thus 



"^0 = T^OJ^fJ' 
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V 



mo. 



\ 



mo,s J 



(2) 



Generators of the u(3) algebra A" = {A° = y2/3/, A*} (where / is an identity matrix, 
A* are Gell-Mann matrices of the su(3) algebra) obeying Tr(A"A^) = 26ab- Due to this 
normalization of this algebra the coupling constants Gg and G^, can be redefined and written 

as Gg = 2/3 Gg, G^ = 2/3 Gi,. 

TABLE 1: Parameter sets of the NJL models. 





NJL (su(2))UJ 


NJL (I)[ . ] 


NJL(II)[ ] 


ENJL [22] 


rriu = rud 


5.5 MeV 


5.5 MeV 


5.5 MeV 


3.61 MeV 







140.7 MeV 


132.9 MeV 


88.0 MeV 


A 


631 MeV 


602.3 MeV 


631.4 MeV 


750.0 MeV 




2.19 


3.67 


3.67 


3.624 


K 





12.36 


9.40 


9.40 


GA^ 











3.842 



The NJL model is non-renormalizable, thus it is not defined until a regularization 
procedure has been specified. This cut-off limits the validity of the model to momenta well 
below the cut-off. In most of our calculations we will adopt the parameters presented in Table 
1. With A, Gg specified above, chiral symmetry is spontaneously broken in vacuum. 

The model contains eight parameters of the standard NJL model (the current mass mo,^, 
mo,s of the light and strange quarks, the coupling constants Gg, determinant coupling K and 
the momentum cut-off A) and additional constant G„ (Gy = XyGg, Xy = 1.06 [12]). In 
the quark massless limit the system has a U{3)l x U{3)r chiral symmetry. The system has 
following global currents: the baryon current 

Jb = l<i^Vq (3) 
and the isospin current which exists only in the asymmetric matter 

J3 = l<iY>^'q , Js = l^Y^'q (4) 

The conserved baryon and isospin charges are given by the relations: 

Q, = ^J d^xq+X'q, z = {0,3,8} 
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which are connected to commuting Cartan algebra. The physical system is defined by the 
thermodynamic potential [23] 

VL = -kTlnTrie-'^^^-'''^^^) (5) 

where H stands for the Hamiltonian. More convenient is to use chemical potentials connected 
to the quark flavor / in the such way that fi^Qi = J2f l^fQf 

/is = - ^/i' (6) 

Quarks and electrons are in /5-equilibrium which can be described as a relation among their 
chemical potentials 

where fiu, jJ^d, li-s and /ig, stand for quarks and lepton chemical potentials, respectively. 
These conditions means that matter is in equilibrium with respect to the weak interactions. 
If the electron Fermi energy is high enough (greater then the muon mass) in the neutron star 
matter muons start to appear as a result of the following reactions 

s ^ u + jjT + 
The neutron chemical potential is 

/in = /i« + 2/id. 

In a pure quark state the star should to be charge neutral. This gives us an additional constraint 
on the chemical potentials 

1 

3"" 3'"" 3' 

where nj (f e u,d, s), Ue are the particle densities of quarks and electrons, respectively. The 
EoS can now be parameterized by only one parameter, namely the dimensionless u quark 
Fermi momentum x {kp^u = M x (M = 939 MeV is the nucleon mass)). 
The mean field approach means that the quantum correlations 

{A- < A >)iB- <B>) = AB-A<B>-B<A> + <A><B>r^O 
may be neglected. It allows us to substitute AB by 

ABr^A<B>+B<A>-<A><B> 
Using this approximation the Lagrange function C may be expressed as 

a=0 

- 2/^E(^/?/) n < ^r^r > +4/i n < m > (8) 



7r^« - - n^s -ne-n^ = 0. (7) 
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D, = d, + -ig.V^X" (9) 



where the co variant derivative is given by 

1 

2' 

Here, the meson fields first appear as non-dynamical variables 

Gs < qX^q >= gsa" (10) 

< qX'^^^q >= g,V^ (11) 

The similar pattern may be extended to the axial mesons 

< qX'^t-f.q >= g^r (12) 
G^ < gA"7^z75g >= g^A"^ (13) 

Mesons masses are defined as 

m„ = gs/\fG~ (14) 

= g^/2\lG~i, (15) 

This is a process of bosonization in which the NJL model produce essentially the m(3) linear 
sigma model as an approximate effective theory for the scalar and pseudoscalar meson sector 
[24]. 

In this paper the variational method based on the Feynman-Bogolubov inequality [25] is 
incorporated (see more details in [ ]) 

VL<VLi = Qo{m^ff)+ < H - Hq >o . (16) 
with the trial Lagrange function described by 

Coimeff) = qih^D^ - meff)q (17) 
and suggested by the mean field form of the Lagrange function (8). The covariant derivative 

D^ = d, + \ig.V^X'' (18) 

is limited to the commuting Cartan subalgebra A* = { A°, A^, }. This approach introduces 
the fermion interaction with homogeneous boson condensates cr", V"^ which together with 
the effective masses mg// will be treated as variational parameters. VLq is the thermodynamic 
potential of the effectively free quasiparticle system 



/n(l + e-/3(v/'='+-e//,/+A'/))) (19) 

We could not forget about energy of quantum fluctuations because it depends on the quark 
effective mass. As fermions give — l/2hujio the vacuum energy, we get 

M rA . 

^o = -;ri E / dkeJ¥+^^f (20) 
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assuming that if rueff = mo the energy of quantum fluctuations may be neglected. The 
effective quark masses entering into the Lagrangian function Lo(me//) of the trial system are 
calculated from the extremum conditions 



dm, 



effj 



(21) 



which give 



or 



where 



ii^eff)fj' = rrieffjSfj' = mc5fj> - 

8 

«=o f'+S 

^eff,f = m,cj -Gs < qjQf >o +2K J| < q^qf >o, (22) 

/V/ 



mfffnNr, k'^dk 



1}. (23) 



exp{/3{^k^ + mlffj + fif)) + 1 
At T = we have only 

/-A 

< ^/g/ >o= -"^e//,/ I dk =■ (24) 

In vacuum we get the constituent quarks with mass 

mc = mo- 2Gs < qq >ov +2K ]J < qf>qf> >o„ (25) 

f'={u,d,s} 

where 

N k'^ 

< QfQf >ov= -^vjT^ / dk = (26) 

^TT^ Jo + ml J 

At minimum the effective free energy has the form 

^e// = ^l\min = ^o("^e//) + Beff (27) 

with the effective bag constant 

Beff = ^Gs < qX^q -4Ji [] < ^/?/ > 

/={n,d,s} 

-^G^ < gT'^A'^g >< q^.X'q > -Bo (28) 

Quarks as effectively free quasiparticles in vacuum with nonvanishing bag 'constant'. The 
constant Bo was chosen in this way to have massive (mv^u,d = 367.61 MeV, = 549.45 
MeV for the NJL (I) parameters set and m^,^u,d = 366.13 MeV, m„,^ = 504.13 MeV for 
enlarged NJL model). However, in high density medium they are less massive (Figure 1) but 

1 /4 

the effective bag constant (Figure 2) grows up to Bjj:j 150 — 180 MeV. The frequently 
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used case with current quarks and bag constant is valid only in very high density limit. This is 
a justification when the quark matter phase is modeling in the context of the MIT bag model 
[9, 27, 28] as a Fermi gas of u, d, and s quarks. In this model the phenomenological bag 
constant By^yr is introduced to mimic QCD interactions. In the original MIT bag model the 
bag constant was constant and the value B = Be = (154.5 MeVY makes the strange matter 
absolutely stable. 

To avoid quantum fluctuations the meson fields may be redefined to produce the 
phenomenological sigma field as 

= Gs{< qX^q >o - < gA'^g >o^,) (29) 
so, the effective quark mass can be rewritten as 

m^ff = m^- g^ipaX" (30) 

Thus, the effective quark mean field (QMF) theory appears. The Lagrange function in the 
mean field approximation may be written in the following form 



(32) 



Unfortunately, the vacuum quantum fluctuations are missed, but meson fields gain the 
dynamical character. Restricted ourself only to u(2) x u(l) subalgebra (a = {0, 1, 2, 3, 8}) 
case we have the simplest version of the QMF theory. Defining new base with r°, a = 
{0,1,2,3,4} as 



for T^ = V2 



/ 1 \ 

r° = 10 t' = 
^ y 

= {1, 2, 3}) the meson fields may be decomposed as follows 




/ \ 




(33) 



and 



Now the new meson fields are denoted by 



UJ 




and fc^ = (with i = {1, 2, 3}), respectively. The simplest u(2) version (cr^ = 0, cj, 
has the Lagrange density function with the following form 



0) 



QMF 



2 



U{a) 



1 

2 

q{iYDf, -mc)q + g^aqq + glSiqr'q 



(34) 
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The field tensors i?" , Vt^y and the covariant derivative are given by 



Rlu = d,K - dub; + gpSa^^b'^b] 



(35) 




(36) 



(37) 



The (5 = limit gives the simplest version of the Quark Mean Field (QMF) model. 

Some years ago Guichon proposed an interesting model concerning the change of the 
nucleon properties in nuclear matter (quark- meson coupling model (QMC)) [30]. The model 
construction mimics the relativistic mean field theory, where the scalar a and the vector meson 
oj fields couple not with nucleons but directly with quarks. The quark mass has to change from 
its bare current mass due to the coupling to the a meson. More recently, Shen and Toki [31] 
have proposed a new version of the QMC model, where the interaction takes place between 
constituent quarks and mesons. They refer the model as the quark mean field model (QMF). 
In this work we shall also investigate the quark matter within the QMF theory motivated by 
parameters coming from the enlarged Nambu-Jona-Lasinio (ENJL). Enlargement of the NJL 
model is based on inclusion of vector mesons while QMF model includes vector mesons at 
the beginning. 

Here the QMF model is a bit generalized by the inclusion of the isovector 5 (ao(980)) 
meson. It splits u and d masses (or proton and neutron masses in the case of RMF approach 
[35, 36]). Both 5* and 6^ mesons may be neglected in the case of symmetric nuclear matter. 
Its role in the asymmetric nuclear matter of the neutron star is significant and is a subject of 
our current interest. 



The QMF model is more flexible. The S'f/(3) symmetry restrict gl = y2/3gs, g^ = g^j, 
gp = gp to gp = gu = gv and mp = m^ = m„. 

The effective quarks mass nipj (or Sf = mpj/M ) dependence on the dimensionless 
Fermi momentum x = kp>a/M is presented on the Figure 1. 

There is not quarks confinement both in the NJL and QMF models. There is no 
mechanism (exept for the NJL solvable model [ ]) to prevent hadrons to decay into free 
constituent quarks. Free constituent quark will produce nearly the same density and pressure 
like free nucleons Sm^ „ ^ ~ M. Without any mechanism of confinement the quark star for 
small densities will have properties very similar to those of neutron stars (the case > 0.65 
in paper [12]) or even white dwarfs. However, this is rather unphysical artefact. It is natural to 
assume that quarks are not allowed to propagate over the distance A ~ "^fT//- this language 
the confinement mechanism introduces infrared cut-off A [ ]. The quarks confinement 
mechanism in the form of the harmonic oscillator potential [ ] may give the nucleon mass 
M = M (a) = M — gNcrC^ + ••• and generate the Relativistic Mean Field (RMF) approach. 

Relativistic Mean Field Theory [32] is very useful in describing nuclear matter and finite 
nuclei. Recent theoretical studies show that the properties of nuclear matter can be described 
nicely in terms of the Relativistic Mean Field Theory. Properties of the neutron star in this 
model were also examined [10, 9, 33, 26, 34]. 





Figure 1. The effective quark masses for different NJL models as a function of the u quark 
dimensionless Fermi momentum x = kp^u/M (M = 939 MeV is the nucleon mass) 
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Figure 2. The effective bag constant as a function of the u quark dimensionless Fermi 
momentum x = kp^u/M (M = 939 M eV is the nucleon mass). 



Its extrapolation to large charge asymmetry is of considerable interest in nuclear 
astrophysics and particulary in constructing neutron star model where extreme conditions 
of isospin is realized. The construction of neutron star model is based on various realistic 
equations of state and results in the general picture of neutron stars interiors. Thus the proper 
form of the equation of state is essential in determining neutron star properties such as the 
mass-range, mass-radius relation or the crust thickness. However, the complete and more 
realistic description of a neutron star requires taking into consideration not only the interior 
region of a neutron star but also the remaining layers, namely the inner and outer crust and the 
surface. The Lagrangian of the RMF theory which tends towards the construction of neutron 
star model contains baryon and meson degrees of freedom and as input quantities coupling 
constants of the mesons and parameters of the potential U (a) which are determined from 
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Figure 3. The binding energy Eb for the quark and nuclear system. The pure quark system 
{usd quarks with the bag constant B) and QMF have nonreahstic behavior for small densities 
as i? 7^ 0. To compare the symmetric nuclear matter (the RMF approach with TM1[39] 
parameterization) is presented. 
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Figure 4. The form of equation of state (EoS) for quark ( NJL(I), enlarged NJL and QMF 
model) matter and nucleon one (TM1[39]). The blue (n^ = 3.94 n^) points the strange s 
quark appearence in the strange star. 



nuclear matter properties. 

Cqmf = \d^od^o - U{a) - ^fi^.fi^'^ + ^M^cu^u'^ + 

i^iiYD^ - M)ij + QnM^ (38) 

Now ^l> describes nucleons and Qn^ = 3(71, Qnw = Sg^ and Qnp = gp,- 

The potential function U (a) may possesses a polynomial form introduced by Boguta 
and Bodmer [37] in order to get a correct value of the compressibility K of nuclear matter at 
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x=kpyM 



Figure 5. The electron and quarks dimensionless Fermi momentum as a function of the u 
quark one {x = kp^^ jM, M = 939 MeV is the nucleon mass). The muon distribution is not 
visible in this scale. 



saturation density 

\J{a) = ^mla' + ^g^a' + ^g^a^ (39) 

Different parameter sets give different forms of the equation of state in the high density 
region above saturation. In this paper the TMl [39] parameter set was exploit. The RMF 
model leads its own fenomenological life but its parameters should be connected to the 
enlaged NIL model. The TM 1 parametrization suggests that x„ = 0.65. 

The p meson plays decisive role in accounting for the asymmetry energy of nuclear 
matter thus its inclusion in a theory of neutron star matter is essential. Also the proton number 
density is determined by this meson. The results for the binding enrgy is presented in Figure 
3. It shows that the nuclean or quark matter in /^-equilibrium has a larger energy per particle 
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than symmetric nuclear matter. For neither parameter set are these matters self-bounded. 
Figure 3 depicts the binding energy for different models. It reproduces the standard results 
for symmetric nuclear matter too. The asymmetric matter is less bound than the symmetric 
one. In this model we are dealing with the electrically neutral neutron or quark star matter 
being in /3-equilibrium. Therefore the imposed constrains, namely the charge neutrality and 
/5-equilibrium, imply the presence of leptons. 



The quark star properties 



To calculate the properties of the quark star we need the energy-momentum tensor. The 
energy-momentum tensor can be calculated taking the quantum statistical average 



where 



T 



2- g^uC. 



(40) 



(41) 



In case of the fermion fields it is more convenient to use the reper field e" defined as follows 



gfiiy = e'^e^riab where rjab is the flat Minkowski space-time matrix. Then 



(42) 



We define the density of energy and pressure by the energy - momentum tensor 

( e = c^p \ 



V 









POO 
P 
OOP/ 



(43) 



where is a unite vector {u^u^ = 1). Both e and P depend on the quarks chemical potential 
H or Fermi momentum Xf. Fermions (quarks and leptons) contribution to the energy and 
pressure are 

" " (44) 



= J2 (^oXsix f,T)+ (^oXL{xf,T) 

f={u,d,s} /={e,/i} 

Pf= E Po<P{xf,T)+ E Po0(x/,T) 

f={u,d,s} f={e,^i} 



(45) 



The fact that effective quarks mass depends rrieffj = 5fM on fermion concentration ( or 
quark chemical potential /i/) now must be included into xi^f, T) and T), 



3 r^/A^ 
x{x,T) = ^ dz 

J A 



Z \ z 



eM{^^\x) + - l^')/r) + 1 



(46) 



exp{U6^{x) + z^ + fi')/T) + l 
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<Pix,T) = -J { , (47) 

} 



where r = {kBT)/M, 



^' = ^/M = ^6^{x)+x^ (48) 

and 

x = k/M (49) 

for each flaver /. Similar to paper [33] we have introduced (48,49) the dimensionless "Fermi" 
momentum even at finite temperature which exactly corresponds to the Fermi momentum at 
zero temperature. To avoid free quark contributions to the equation of state coming from 
small densities the infrared cut-off \ = 5 [40] was introduced. The case A = with the NJL 
(I) parameters set nicely reproduce the result od paper [12]. 

The parametric dependence on // (or Xf) defines the equation of state. The various 
equations of state for different parameters sets are presented in Figure 4. The binding energies 

Eb = p/uB - Mc' 

for the bulk nuclear {ub is the baryon number density, ub = {rip + n„)) and quark matter 
(ub = {riu + nd + Us)/ 3) are presented in Figure 3. 

The metric is static, spherically symmetric and asymptotically flat 

/ e^(^) \ 

-e^(^) 

-r2 

V -r^sin^e J 

(where z/(r) and A(r) are functions of a radius r ). Einstein equations (in isotropic case) 
leads to the standard Tolman-Oppenheimer-Volkoff equations [-fi]. The equations describing 
masses and radii of quark star are determined by the proper form of the equation of state. The 
obtained form of the equation of state is the base for calculating macroscopic properties of the 
star. In order to construct the mass-radius relation for given form of the equation of state the 
OTV equations have to be solved 

^ = Anr^pir) (52) 
dr 

The continuity condition for the energy-momentum tensor T'^'^ = defines the connection 
between gravitational potential ^{r) (50) and pressure and density profiles P(r) and p(r) 
du{r) 2 dP{r) 



(50) 



dr P{r) + c^p{r) dr 



(53) 
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The equation (52) determines the function A(r) 

r 

Having solved the OTV equation the pressure P(r), mass m(r) and density p(r) profile is 
obtained. To obtain the total radius R of the star the fulfillment of the condition P{R) = is 
necessary. Introducing the dimensionless variable ^, which is connected with the star radius 
r by the relation r = enables to define the functions P(r), p(r) and m(r) 

p(r)=poxWO) (54) 

P(r) = PofixiO) (55) 

m(r) = Mqu{0 (56) 

by ^. Dimensionless functions defined as 

"0 = ^ , /5o = 3—, M, = -vrpoa (57) 

Poa Mq 3 

are needed to achieve the OTV equation of the following form 

— = -ao(x(a^(0) + ^{Ai))) ^^^^ 

^ = /5ox(x(0)e' (59) 

with being the gravitational radius. The equations (58,59) are easy integrated numerically. 
These are equations for dimensionless mass m(^) = m(r)/M0 up to radius dimensionless ^ 
and the quark u dimensionless Fermi momentum x = kp^u/M. Knowing the variable x all star 
properties can be calculated. Quarks and electron dimensionless Fermi momenta dependences 
on X is presented in Figure 5. 

Both nuclear and quark matter being in /^-equilibrium are not bound (Figure 3). For 
quark matter at moderate densities is bound due to the presence the bag constant B which acts 
as a negative pressure P = —B + ... (Figure 4).Higher density matter is bound only gravity. 
The gravitational binding energy of the star is defined as 

E^^g = {Mp - m{R))c^ (60) 

where 

M, = 47r [""drr^l - '^-^^^Ylp[r) (61) 

is the proper star mass. 

Strange quarks star in the NJL model are rather small in comparison to the neutron ones. 
To compare these strange stars to neutron stars models obtained in the RMF approach the 
mass - radius relations are also presented in Figure 6. ( the black solid line (TMl) for pure 
npe nuclear matter and the star with a crust (TMl + Bonn + Negele + Vautherin [42] - dotted 
line). The smaller sizes of quarks star is due to the fact that the pressure (Figure 4) reaches 
zero for high densities {rf§ ~ 0.26 fm~^ = 1.75 ra^) then to the symmetric nuclear matter 
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{n% ~ 0.15 fm~^). On the surface such a star has higher density than the saturated nuclear 
matter. This makes the star smaller and denser. The same situation exists in the QMF approach 
when the bag constant is included. In the NIL model the effective bag constant has dynamical 
origin (eq. 28). Its main contributions come from quantum vacuum fluctuations and scalar 
mesons. The strange quark s appears only for densities high enough {ub > n% = 3.94 n^). 
For smaller densities the quark star is built only from u and d quarks. It is important to stress 
that even then the quantum vacuum fluctuation come from all flavors including s quark and 
all antiparticles. There are no s quarks for small densities in the strange quark star only its 
quantum fluctuations. In the QMF approach there is no quantum fluctuations at all. This is a 
significant difference between the NIL and QMF approach. The another one is that the bag 
constant in the QMF model must be added by hand. This make the QMF approach unreliable 
for smaller densities (Figure 4). 

The gravitational binding energy of a strange quark star and QMF approach) and neutron 
(RMF) star is presented in Figure 10. The arrow shows a possible transition from unstable 
neutron star to the strange one with conservation of the baryon number Mb = Modot, {Mb = 
rrinC^NB). 

In the enlarged NIL model [22] vector mesons are included. There contribution to the 
effective bag constant are positive (Figure 2). This makes that a strange quark star possesses 
a bit bigger radius and mass then in the NJL (I) model. 

However, a maximum stable strange quark star is obtained for p2 = 3.1 10^^ g/cm^ 
and has the following parameters M = 1.61 Mq and R = 8.74 km. The baryon number 
is this star is the same as in the case of pure neutron star with Mb = 2.126 M©. Below 
the density ps = 3.94 po (Po = 2.5 10^^ g/cm^), there are no strange quarks and quark stars. 
Stable stars are those with ^ > [38] (Figure 7). A gravitational binding energy for a 
strange quark is lower then a neutron one for p > 1.6 10^^ g/cm^. The M(p) dependence for 
the quark star is presented in Figure 7. For the quark star with the maximal central density 
Pc = 3.11 10^^ g/cm^ the star profile is presented on the Figure 8. Quark and electron mass 
distributions inside the star are presented in Figure 8. The quark partial fractions defined as 

X = ^/ _ ^ 

^ {nu + nd + Us) ub ' 

where / = (u,d,s) are presented in Figure 5 . 
Conclusion 

The properties of the strange quark star in the bag model with B = Be and the current quark 
masses (m„ = = 0, = 150 MeV) are presented in [9]. The star model base on the 
QMC is very similar (the dotted green line in Figure 6) and close to the one for the ENJL 
model parametrized by the TMl parameter set (x^ = 0.65, the solid violet curve 6). The 
properties of strange stars with quark masses changing continuously from constituent quark 
mass to the small current (see Figure 1) are presented in [43]. All these stars are more compact 
then neutron stars (see the Figure 6) and are similar to those of the NJL (I) model. 
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Figure 6. The mass - radius M{R) dependence for the quark star(NJL(I), ENJL (x„ = 1.06, 
the solid brown line and Xy = 0.65, the solid violet curve), QMF model and pure uds matter 
(QMC)). To compare this relation to the neutron star in the RMF approach it is presented as the 
black line (TMl - solid line) for pure npe nuclear matter and (TMl+Bonn+Negele+Vautherin) 
- dotted line) with crust. 



In this paper the enlarged NJL model is used to construct the EoS and properties of the 
strange quark star. The stable strange quark star exists from minimal central density till the 
maximal one p2 = 3.1 10^^ g/cm^ which gives the following star parameters M = 1.61 Mq 
and R = 8.74 km. Its baryon number is the same as for the pure neutron star with 
Mb = 2.126 Mq. The very similar strange star but less compact is obtained in the solvable 
NJL model [13]. The gravitational binding energy for a strange quark star is lower then the 
neutron one for densities p > 1.6 10^^ g/cm^. The conversion of an unstable neutron star 
into a strange star is an exciting subject which may help explain the gamma-ray bursts enigma 
[44]. 
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Figure 7. The dependence for the quark star M{pc) mass on the its central density pc- 



Similarly to the QMF model the enlarged NJL one includes also the coupling to vector 
mesons. This is crucial for the quark star properties. Also the quark s mass is important. The 
mass of a quark s is also relevant because its smaller mass causes that a strange quarks appear 
for lower densities. Nonzero strangeness of the matter gives as a result a strange star. It is 
fascinating that the neutron and strange star properties are strictly connected to inner structure 
of nuclei and nucleons. 
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Figure 8. The profile for effective quark Fermi momentum of the maximal quark star (ENJL 
model) with the central density pc ~ 3.11 10^^ g/cm?. 
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